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Summary 


It  is  known  that  the  type  of  the  system  of  partial  differential 
equations  governing  finite  elastostatics  can  change  from  elliptic  to  non- 
elliptic  at  sufficiently  large  deformations  for  certain  materials.  This 
introduces  the  possibility  that  the  elastostatic  field  may  exhibit  certain 
discontinuities.  Some  aspects  of  the  general  theory  associated  with  these 
issues  were  examined  in  a recent  series  of  studies  by  Knowles  and 
Sternberg.  In  this  paper  we  illustrate  the  occurrence  of  elastostatic 
fields  with  discontinuous  deformation  gradients  in  a physical  problem. 

The  body  is  assumed  to  be  composed  of  a material  which  belongs  to  a 
particular  class  of  isotropic,  incompressible,  elastic  materials  which 
allow  for  a loss  of  ellipticity.  It  is  shown  that  no  solution  which  is  smooth 
in  the  classical  sense  exists  to  this  problem  for  certain  ranges  of  the 
applied  loading.  Next,  we  admit  solutions  involving  elastostatic  shocks 
into  the  discussion  and  find  that  the  problem  may  then  be  solved  completely. 
When  this  is  done,  however,  there  results  a lack  of  uniqueness  of  solutions 
to  the  boundary-value  problem.  In  order  to  resolve  this  non- uniqueness, 
the  dissipativity  and  stability  of  the  solutions  are  investigated. 


Hie  results  communicated  in  this  paper  were  obtained  in  the  course  of  an 
investigation  supported  by  Contract  N00014-75-C-0196  with  the  Office  of 
Naval  Research  in  Washington  D.  C. 


1 . 1 Introduction 


In  the  course  of  investigating  some  crack  problems  [l],  [2],  [3], 
Knowles  and  Sternberg  encountered  certain  difficulties  which  suggested 
that  the  problem  may  not  admit  a classically  smooth  solution.  In  order 
to  clarify  this  situation,  a series  of  preliminary  studies  were  undertaken 
(References  [4]  - [ 7 3 > in  which  these  authors  looked  at  the  question  of 

the  change  of  type  of  the  displacement  equations  of  equilibrium  from  elliptic 
to  non-elliptic,  and  the  related  issue  of  the  existence  of  solutions  possessing 
certain  discontinuities  — referred  to  as  elastostatic  shocks.  The  presence 
of  such  elastostatic  shocks  is  found  to  affect  the  energy  balance  of  the 
field.  This  led  Knowles  and  Sternberg  [6],  [7]  to  propose  a notion  of 
dissipativity  associated  with  such  fields . Subsequently,  Abeyaratne  [8] 
examined  the  corresponding  issues  in  the  case  of  incompressible  materials. 

In  order  to  illustrate  the  occurrence  of  elastostatic  shocks  in  a 
boundary -value  problem,  we  consider  a problem  in  finite  plane  strain  for 
a hollow  circular  cylinder.  Specifically,  we  examine  the  case  in  which 
the  outer  surface  of  the  cylinder  is  held  fixed  while  the  inner  surface  is 
twisted  circumferentially.  The  cylinder  is  presumed  to  be  composed  of  a 
homogeneous,  isotropic,  incompressible  elastic  material.  Although  this 
problem  has  been  considered  before*,  our  interest  centers  on  those  materials 
whose  strain  energy  density  permits  a failure  of  ellipticity  of  the  displacement 
equations  of  equilibrium  at  sufficiently  severe  deformations. 

We  demonstrate  for  our  choice  of  material  that,  while  for  both 
sufficiently  large  and  small  values  of  the  prescribed  twist  the  problem 


*See  Rivlin  [9]  as  well  as  Green  and  Zerna  [10]  page  95.  The  problem 
has  been  reconsidered  more  recently  by  Ogden,  Chadwick  and  Haddon  [ 11  ]. 
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admits  a unique  smooth  solution,  there  are  certain  intermediate  ranges  of 
the  prescribed  twist  at  which  no  classically  smooth  solutions  exist.  We 

% 

then  show  that  there  are  however,  an  infinite  number  of  weak  solutions 
• involving  elastostatic  shocks  in  these  ranges  of  the  applied  twist. 

We  then  consider  the  quasi-static  problem  in  which  the  prescrib- 
ed twist  is  gradually  changed  in  time,  and  explore  the  consequences  of  the 
dissipation  inequality.  It  turns  out  that  enforcing  this  inequality  fails  to 
single  out  a unique  weak  solution. 

In  an  attempt  to  clarify  this  issue  of  non- uniqueness,  we  examine 
the  stability  of  the  various  equilibrium  solutions  against  purely  circum- 
ferential perturbations.  It  is  found  that  the  classical  energy  criterion 
for  stability,  without  reference  to  the  dissipation  inequality,  picks  out 
t a unique  solution  to  the  boundary -value  problem  at  every  value  of  the 

prescribed  twist.  In  the  discussion  of  the  various  issues  outlined  above, 
we  restrict  attention  to  configurations  involving  not  more  than  one  elasto- 
static shock.  As  a consequence  of  the  stability  criterion,  we  find  that  an 
equilibrium  solution  involving  more  than  one  shock  is,  in  fact,  unstable. 

Ericksen  [ 12]  has  discussed  the  equilibrium  of  a bar  composed 
of  a material  whose  stress  response  in  uniaxial  tension  is  qualitatively 
similar  to  the  shear  stress  response  in  simple  shear  of  the  class  of 
materials  considered  here.  There  is  a striking  similarity  between  his 
results  and  ours;  in  fact,  certain  aspects  of  our  study  of  the  stability 
of  weak  solutions  were  suggested  by  the  arguments  in  [ 12], 

In  Section  2 we  set  up  the  classical  problem  governing  the  twist- 
i ing  of  a hollow  cylinder  composed  of  an  arbitrary  homogeneous,  incompres- 

sible, isotropic,  elastic  solid.  We  then  discuss  the  particular  class  of 
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materials  with  which  we  will  be  concerned.  In  Section  3 we  determine 
the  solutions  of  this  problem  and  construct  the  associated  torque 'twist 

% 

curves.  For  sufficiently  small  and  large  values  of  the  prescribed  twist, 
we  have  a unique  smooth  solution  at  which  the  displacement  equations  of 
equilibrium  are  elliptic.  Depending  on  the  details  of  the  geometry  and 
constitutive  law,  it  is  also  possible  to  have  a unique,  non-elliptic,  smooth 
solution  at  certain  — but  not  all  — values  of  the  twist  in  the  intermediate 
range.  In  all  cases  there  are  ranges  of  values  of  the  prescribed  twist 
for  which  we  find  no  solution.  We  then  prove  that,  in  fact,  no  smooth 
solutions  exist  in  these  ranges  of  the  prescribed  twist. 

We  next  set  up  and  solve,  in  Section  4 the  problem  in  its  weak 
formulation.  We  now  find  a solution  corresponding  to  every  value  of  the 
t prescribed  twist,  but  unfortunately  there  are  many  solutions  corresponding 

to  certain  twist  values. 

% 

In  Section  5 we  make  use  of  the  dissipation  inequality  in  an  unsuc- 
cessful attempt  to  extract  a unique  solution  from  among  the  many  solutions 
to  the  boundary  value  problem.  Finally,  in  Section  6 we  examine  the 
stability  of  each  of  the  available  solutions  against  purely  circumferential 
perturbations.  We  find  that  at  every  value  of  the  prescribed  twist  there 
is  precisely  one  stable  solution  to  the  boundary -value  problem  in  its  weak 
formulation.  For  sufficiently  small  and  large  values  of  the  applied  twist. 


/ 


this  unique  stable  solution  is  smooth  and  elliptic.  For  all  intermediate 
values,  the  unique  stable  solution  involves  an  elastostatic  shock  and  is 
elliptic. 


-5- 


2.  1 Formulation  of  Problem 

Suppose  that  the  open  region  ft  occupied  by  the  interior  of  a 
body  in  its  undeformed  configuration  is  a hollow  right  circular  cylinder 
of  internal  and  external  radii  a and  b , respectively.  Let  II  be  the 
open  middle  cross-section  of  the  cylinder  ft  , and  let  O be  the  center 
of  the  annular  region  II  . 

Suppose  the  inner  surface  of  the  cylinder  is  rotated  circumfer- 
entially through  an  angle  , while  the  outer  lateral  surface  is  held 
fixed.  We  assume  that  the  resulting  deformation  maps  the  point  with 
cylindrical  coordinates  (r,  8,  z)  in  the  undeformed  configuration  onto 
the  point  with  cylindrical  coordinates  (p,  if,  5),  where 

p = p(r,  6,  z)  = r , t = *(r,  8,  z)  = Q + 0(r)  , S =§(r,  8,z)  = z . (2.1) 

This  describes  a plane  deformation  in  which  each  particle  moves  cir- 
cumferentially through  an  angle  0(r).  Suitable  tractions  are  presumed 
to  be  applied  on  the  ends  of  the  cylinder  so  as  to  maintain  such  a state 
of  plane  strain. 

The  deformation  (2.1)  may  be  equivalently  expressed  as  follows. 
Let  X be  a fixed  rectangular  cartesian  coordinate  frame  with  its  ori- 
gin at  O and  ortho  normal  base  vectors  e.  , e,  , e-  , such  that  e,  and 
£2  are  in  the  plane  of  II  and  £3  is  normal  to  II  . If  £ is  the  position 
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vector  after  deformation  of  the  particle  which  was  located  at  x in 
the  undeformed  configuration,  we  can  write  (2.  1)  as 


where 


y-pXjCOs  0(r)  - x2sin  ^(r)  , 
y 2=  x2cos  0{r) + XjSin  0(r)  , > 

y3=X3  ' 

J 


r 


(2.2) 


(2.3) 


Here  y.  and  x.  , i=l,2, 3,  are  the  components  of  the  vectors  y and 
x in  the  frame  X . We  will  temporarily  assume  that  the  local  angle 
of  twist  0(r)  is  twice  continuously  differentiable  on  (a,  b). 

It  is  convenient  to  express  the  field  quantities  at  any  point 
(r,9,z)  in  terms  of  components  in  the  rectangular  cartesian  coordinate 
frame  X*  which  is  obtained  by  rotating  the  frame  X through  an  angle 
f(r,  0,  z)  about  the  e.-axis.  The  matrix  of  components  of  the  deforma- 
tion  gradient  tensor  F=Vy  in  the  frame  X',  F , is  easily  computed 
from  (2.2),  (2.3)  and  the  change -of -frame  formula  for  tensors  to  be 


X' 
F*  = 


cos  0(r) 


■sin0(r) 


sin  0(r)  + r0'(r)cos  0(r)  cos  0(r)  - r^rjsin  0(r)  0 


(2.4) 


X 1 

Note  that  the  matrix  F may  be  decomposed  as  follows; 
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1 0 0 cos  0(r)  -sin^(r)  0 

v ' 

F = r #'(r)  1 0 sin0(r)  cos0(r)  0 , 


(2.5) 


L 0 0 ‘J  l 0 0 

which  implies  that  locally  the  deformation  (2.  1)  is  composed  of  a rigid 
rotation  through  an  angle  0 about  the  - axi s followed  by  a simple 
shear  parallel  to  the  circumferential  direction  with  an  amount  of  shear 
r0'(r).  Set 


k(r)  = r0'(r)  , 


(2.6) 


so  that  k(r)  is  the  local  amount  of  shear. 

Equation  (2.4)  indicates  that  detF  = 1 , ao  that  the  deformation 

(2.1)  is  locally  volume  preserving.  From  (2.4)  and  (2.6)  we  have 

T 

the  components  of  the  left  Cauchy-Green  tensor  G = FF  : 

f 1 Mr)  ol 


GX=  k(r)  l + k2(r)  0 


(2.7) 


The  principal  invariants  of  G are  found  from  (2.7)  to  be 


I1=  trG  = 3 + k*(r)  , 

I2=l{«trS>2-«S2}=5  + l‘ZW  . 


I,  = det  G = 1 . 
3 ~ 


(2.8) 


Suppose  that  the  body  is  composed  of  a homogeneous,  isotropic. 
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incompressible,  elastic  solid  which  possesses  an  elastic  potential 
W = W(I,,I,).  W represents  the  strain  energy  density  per  unit  unde- 
formed  volume.  The  constitutive  law  for  the  Cauchy  stress  tensor 
t is  then 


,/8Wxr  awv 

i=2v^  + I p-i 


(2.9)1 


where  p(y)  is  a pressure  field  arising  because  of  the  constraint  of  in- 
compressibility.  We  suppose  for  the  moment  that  p(y)  is  continuously 
differentiable  on  ft  . Using  (2.7),  (2.8)  and  (2.9)  we  find  that  the 
stresses  induced  by  the  deformation  (2,  1)  are  given  by 


T^,=  2W/(I)-q  , 

T?2/=T21=2kW'(I)  ' 

T*2'=2(l+k2)W'(I)-q  , 

X'  X'  X'  X' 

T13  31  23  32~°  ' 

* 

T^C,/=2W'(I)  + 2k2  |^(I+ 1,  1+  1)  - q , 


f (2.10) 


33 


where  we  have  set 


I = 2 + k2(r)  , W(I)  = W(I+1,  1+1) 


for  1*2 


(2.11) 


and 


*See  Truesdell  and  Noll  [ 13  ] , page  319. 
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* 

, 8W 
q = P-2^ 


(2.12) 


Since  the  pressure  p might  depend  on  the  coordinates  t and  5,  in 
addition  to  p , it  follows  that  _ne  Cauchy  stress  tensor  jt  might  depend 
on  all  three  of  p , f , and  5 . 

It  is  suggestive  to  introduce  the  notation 


X' 

T = T,  , , 

PP  11  ' 


X' 

22  ' 


(2.  13) 


X ' 

T . = T . = T , _ . 

Pt  *P  12 


The  equilibrium  equation  in  the  axial  direction  is  easily  shown 
to  be  satisfied  if  and  only  if  q does  not  depend  on  y^  (and  hence  5). 
It  follows  from  (2.  10),  that  the  Cauchy  stress  tensor  jr  is  also  inde- 
pendent of  the  axial  coordinate  5 . The  remaining  two  equilibrium 
equations  now  take  the  form 


8t  , 8t  . 

pp  + A _£&  + A(T  . T i = o 

8p  P 9+  P PP  ' 


or  . , 8t , . _ 

^i  + A Ji±  + 2T  0 . 

8p  p 8t  p p* 


(2.14) 


(2.15) 


From  (2.1),  (2.6),  (2.10),  (2.11)  and  (2.13)  we  see  that 
is  independent  of  the  coordinate  if  , whence  (2.  14)  and  (2.  15) 
specialize  to 


‘Body  forces  are  presumed  to  be  absent. 


: • * i 3 


-10- 


8t  t - T 

“P  P ' 


(2. 16) 


(2.  17) 


Integration  of  (2.  17)  with  respect  to  i)r  leads  to 


TM=-{p!?(p2Tp»)}*+<:<l)l  • 


(2.  18) 


where  c(p)  is  a "constant"  of  integration  depending  on  p alone.  It  is 
apparent  from  (2.18)  that  t is  single- valued  only  if 


(2.19) 


It  now  follows  from  (2.18)  and  (2.19)  that  t , and  hence  q,  p and 
Tpp  as  we^»  are  independent  of  the  angular  coordinate  i|r  . Using 
(2.10)  and  (2.13)  in  (2.16)  and  (2.19)  we  obtain  the  governing 
system  of  ordinary  differential  equations  for  0(p)  and  q(p): 


d_ 

dp 


{p30,(p)W'(2  + pV2(p))}  = O , 


(2.20) 


2|^W/(2  + pV2(p))-2p0'2(p)W/(2  + pV2(p))  = -|a  . 


(2.21) 


On  integrating  (2.20)  with  respect  to  p we  find  that 


p30/(p)W/(2  + pV2(p))  = -^  » 


(2.22) 


where  T is  a constant  of  integration.  Likewise,  integration  of  (2.21) 


mm 
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with  respect  to  p and  making  use  of  (2.22)  gives 

P 

q(p)  = 2W/(2  + pV2(p))+^  J ^-dS+qj,  , (2.23) 

a 5 

where  is  a constant. 

It  is  convenient  to  define  the  scalar  valued  function  f by 

f(k)  = 2kW'(2  + k2)  for  - oo<k<oo  . (2.24) 

It  is  readily  seen  that,  if  an  incompressible,  isotropic,  elastic  solid  is 
subjected  to  a simple  shear  deformation,  the  shear  stress  corresponding 
to  an  amount  of  shear  k is  f(k)  . Accordingly  the  function  f may 
be  interpreted  as  the  shear-stress  response  function  in  simple  shear. 
Equation  (2.22)  can  now  be  written  as 

f(p0'(p))  = - — on  (a,  b)  , (2.25) 

2ttP 

which,  together  with  the  boundary  conditions 

0(a)  =0Q  , (2.26) 

0(b)  = 0 , (2.27) 

constitutes  the  boundary  value  problem  for  0(p).  We  wish  to  find  a 
function  0(p),  continuous  on  [a,  b]  and  twice  continuously  differentiable 
on  (a,b)  , and  a real  number  T such  that  (2.25)  - (2.27)  hold.  We 
will  refer  to  such  a solution  as  a smooth  solution.  Note  that  once  0(p) 
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has  been  so  determined,  (2.23)  gives  q(p)  directly. 

Finally,  note  from  (2.6),  (2.10),  (2.13)  and  (2.22)  that 
2 

Tp^(p)  = - T/2irp  so  that  T is  the  torque  per  unit  axial  length  of  the 
cylinder  acting  on  the  inner  surface,measured  positive  in  the  counter- 
clockwise sense. 

2.  2 A Particular  Class  of  Constitutive  Laws 

We  now  describe  the  particular  class  of  homogeneous,  isotropic, 
incompressible,  elastic  materials  to  which  we  will  restrict  attention 
in  this  study.  It  is  adequate  for  our  purposes  to  specify  the  response 
of  the  material  in  simple  shear  alone.  Observe  from  (2.24)  that  W(I) 
is  completely  determined  by  the  function  f . Consequently,  one  can  show 
that  the  response  in  simple  shear  determines  completely  the  in-plane 
response  in  all  plane  deformations  for  such  materials. 

Equation  (2.24)  implies  that  f is  an  odd  function,  i.e., 

f(k)  = -f(-k)  for  - oo<  k<  ao  . (2.28) 


We  assume  that 

(i)  f is  continuously  differentiable  on  (~ao,oo), 

(ii)  f is  positive  on  (0,  oo)  , whence  it  follows  from  (2.28)  that 


kf(k)>0  for  k^O  , 


(2.29) 


(iii)  there  exist  real  numbers  kj  and  k^  (0<kj<k2<oo)  such 


that 

fUjHfU^O  , 

f'(k)>0  for  0sk<kj  , k2<k<OL  , 
f'(k)<0  for  kj<k<k^  , 


>v 


> (2.30) 
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(iv)  f(k)-*oo  as  k -»oo  . 

Since  f is  an  odd  function  it  now  follows  that  f'(-k j)  = f'(-k2)  = 0 , 
f/(k)>0  for  -ao<k<-k2,  -kj<ksO,  f'(k)<0  on  -k2<k<-kj  and 
f(k)-»-co  as  k -*  - oo  . Therefore,  the  function  f(k)  has  local  maxima 
at  k = kj  , -k2  and  local  minima  at  k = k2  , - k ^ and  is  monotone  in 
between.  A graph  of  such  a function  f is  shown  in  Fig.  1,  where  we 
have  set 


f (k - ) = t , f(k,)  = t 

1 max  ' '2'  min 


(2.31) 


Note  that  necessarily 


U(k)|sT 


max 


t . s | f (k)  | £ T 
min  11  max 


|f(k)|i 


min 


for  |k|skj  , 
for  kjS  | k | £ k2  , 
for  | k | i k2  . 


> 


(2.32) 


An  isotropic,  incompressible,  elastic  solid  subjected  to  plane 
deformations  conforms  to  the  in-plane  Baker- Ericksen  inequality  if  and 
only  if  W'(I)>0  for  I>2  . By  virtue  of  (2.24)  this  is  equivalent  to 
kf(k)>0  for  k^O  . Because  of  (2.29),  the  class  of  materials  under 
consideration  satisfies  this  condition. 

Moreover,  we  have  from  Section  3.2  of  [8]  that  in  any  plane 
deformation,  the  plane  strain  displacement  equations  of  equilibrium 
are  elliptic  at  a point  if  and  only  if  the  associated  local  amount  of  shear  * 
is  less  than  kj  or  greater  than  k2  . In  the  context  of  the  problem 
considered  here,  we  have  from  (2.6)  and  (2.  11)  that  Jjjg. 


T 


See  Section  2.2  of  [8]. 
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displacement  equations  of  equilibrium  are  elliptic,  at  a solution  cor 


responding  to  0(r)  and  a point,  if  and  only  if  |r0'(r)|  is  less  than 


It  is  clear  that  for  the  particular  class  of  materials  just  de 


scribed,  f has  no  single-valued  inverse.  The  restrictions  of  f to 


certain  subintervals  of  (-00,00),  on  the  other  hand,  do  have  unique 


inverses.  Let  F,  , F,  and  F,  be  the  functions  defined  by 


By  virtue  of  (2.32)  , it  follows  that  Fj  , F^  and  F^  are  de- 
fined  on  l-Tmax.  *miJ  . [-Tm„  . u ITmin  • 'mJ  *nd 

(-00  , -Tminl  U [ Tm^n  » 00)  respectively,  and  that  they  are  continuously 
differentiable  on  the  corresponding  open  intervals. 

The  following  properties  of  the  inverses  F.  (i  = 1,2,3)  can  be 
easily  verified;  we  list  them  here  for  subsequent  reference. 


min 
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F '(t)>  0 for 

1 T 1 < T , 

1 1 max  ' 

> 

F2(t)<0  for 

t . < It  I <t  , 

min  1 1 max  ' 

► 

(2.35) 

F^(t)>0  for 

Tmin<lTl<”  - 

> 

> 

F.(t  ) 

1 max 

= F2(Tmaxl=kl  ’ 

► 

(2.36) 

F,(T  ): 

2'  max' 

= F,(t  ) = k,  , 

3 max  2 ’ 

J 

1 Fj(t)  | < kj 

for  1 T 1 s T , 

1 1 max  ' 

kjS  |F2(T)|sk2 

,or  TminS  |,|STm.x  ' 

(2.37) 

|F3(x)|ik3 

tor  lTl*Tmi„  • 

y 

F3(|t|)>F2(|t|)>F1(|t|)  for  Tmin<  I T I < Tmax  • <2-38) 

F.(-r)  = -Fi(T)  , 
tF1(t)>0,(t^0)  , 

|f.(t)|  = f4(  It  I ) , 

Fj(t)  -»iCD 


A 


i = 1,  2,  3 and 
for  t in  the 
' appropriate 
interval 


as  t-»±qd 


(2.39) 

(2.40) 

(2.41) 


(2.42) 
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3.  1 Smooth  Solutions 

We  now  return  to  the  task  of  solving  (2.25)  - (2.27)  for  the 

special  class  of  materials  described  in  the  previous  section.  To  this 

end,  we  first  establish  the  following  preliminary  result. 

Lemma:  There  does  not  exist  a solution  0(r)  in 
2 

the  class  C (a,b)  to  the  differential  equation 

f(r*'(r))  = - T/2*r2  (3.1) 

where  f is  a continuously  differentiable  function 
conforming  to  (2.  28)  - (2.30)  and  T is  a con- 
stant, such  that  at  some  radius  s , a<s<b  , 

s^/(s)  = ±k1  or  ±k2  . (3.2) 

Proof:  Suppose  that  there  is  such  a solution  0(r).  Differentiation  of  (3.  1) 
with  respect  to  r and  setting  r = s leads  to 

f/(s0,(s)){0#(s)  + s0'(s)}  = T/Trs3  , 

which  because  of  (2.28),  (2.30)  and  (3.2),  yields 

T=0  . (3.3) 

The  differential  equation  (3.1)  now  reads  f(r$'(r))  = 0 , 


...  -i 

' 
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which  because  of  (2.28)  and  (2.29)  implies  that 

0'(r)  = O on  (a,  b)  . (3.4) 

Equation  (3.4)  , however,  contradicts  the  assumption  (3.2)  . This 
establishes  the  lemma. 

Now  suppose  that  the  prescribed  twist  is  a number  in  the 
interval 


(3.5) 


On  using  (2.26)  , (2.27)  , (2.35)  and  (2.37)  in  (3.5)  we  have 


whence 


(3.6) 


J 


50'(5)  + k, 

r -d5  > 0 


? W'«)-k. 

r -d5<0  . 


a a 

From  the  preceding  lemma  we  know  that  r0'(r)^±kj  on  (a,  b), 

(3.7)  implies 


(3.7) 


so  that 


k1>r0/(r)>- kj  on  (a,  b)  , 


(3.8) 


since  the  integrands  in  (3.7)  are  continuous  on  (a,b).  Therefore,  if 
is  a number  such  that  (3.5)  holds,  then  necessarily  the  solution  to 
(2.  25)  - (2.  27)  must  satisfy  (3.8)  But  because  of  (2.33)  and 
(2.34),  we  see  that,  (2.25)  and  (3.8)  hold  if  and  only  if 
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r*'(r)  = F.(-  -2~z)  . (3-9) 

1V  2irr  ' 

Integrating  this  and  using  the  boundary  condition  (2.27)  together  with 
(2.39)  leads  to 


(3.  10) 


On  enforcing  the  boundary  condition  (2.26)  , we  have  from  (3.  10)  that 


(3. 11) 


Finally,  we  verify  that  (3.11)  determines  a unique  number  T 
for  every  given  number  in  the  interval  defined  by  (3.  5)  . To  this 
end,  define  the  function  by 


for  It  I £ 2 it  a t 


max 


(3.  12) 


On  differentiating  (3.12)  with  respect  to  T and  making  use  of  (2.35), 


we  find  that 


# j(T)>0  for  |T|*2Tra2Tmax  , 


(3.13) 


whence  is  mono  tonic  ally  increasing  on  [-2Tra2Tmax,  2ira  TmajJ  * 

Thus,  if  0O  is  a number  such  that 

^<-2"a2Tmax)S:  V’l^^max*  ’ 


(3.  14) 
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then,  0q=#j(T)  defines  a unique  number  T . Note  that  (3.  14) 
because  of  (3.12),  is  identical  to  (3.5). 

Therefore,  we  conclude  that,  if  the  prescribed  twist  0^  is  in 
the  interval  defined  by  (3.5)  , Equation  (3.11)  determines  a unique 
real  number  T , which  together  with  (3.  10)  gives  the  corresponding 
unique  smooth  solution  to  (2.25)  - (2.27). 

In  an  entirely  analogous  manner,  we  can  show  that,  if  the  pre- 
scribed twist  satisfies 


then,  the  relation 


»o  = J 

a 2tt5 


(3.16) 


determines  a unique  real  number  T , which  together  with 


*<r>  = f?ri(T7l)d5 


Zni 


(3. 17) 


is  the  corresponding  unique  smooth  solution  to  (2.  25)  - (2.27). 
Similarly,  if  the  prescribed  twist  0^  obeys 


(3.  18) 


the  relation 
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*°=HF3^)d! 

determines  a unique  real  number  T , which  together  with 


b 


(3.  19) 


(3.  20) 


is  the  corresponding  unique  smooth  solution  to  (2.  25)  - (2.  27). 

We  will  refer  to  (3.  10)  , (3.  17)  and  (3. 20)  as  (smooth)  Solution  1, 
Solution  2 and  Solution  3 respectively.  Equations  (3.  11)  , (3.  16) 
and  (3.19)  are  the  corresponding  torque-twist  relations.  One  sees 
readily  from  (3.  10)  , (3.  17)  , (3.  20)  and  the  discussion  of  ellipticity 
in  Section  2.2  , that  the  displacement  equations  of  equilibrium  are 
elliptic  everywhere  in  n at  Solution  1 and  Solution  3,  and  that  they 
are  non-elliptic  at  Solution  2. 

Because  of  (2.35) one  has 


Jt 


?*2 


d?  , 


(3.21) 


if  and  only  if 


. 2 2 
b t . s a t 

mm  max 


(3.22) 


Accordingly,  it  is  only  when  (3.22)  holds  that  there  are  values  of  <6^ 
in  the  interval  (3.  15)  , and  consequently  that  Solution  2 exists.  In  this 
paper,  we  will  consider  in  detail  the  case  when  the  dimensions  of  the 
tube  and  the  constitutive  law  of  the  material  are  such  that 


2 T 

a ^ mm 

72>  t ' 

b max 


(3.23) 


and  trace  the  other  cases  (which  are  in  fact  less  complicated)  through 

footnotes.  The  end  result  turns  out  to  be  the  same  in  all  cases.  For  a 

given  material,  one  could  view  (3.23)  as  requiring  the  thickness  of  the 

tube  to  be  sufficiently  small.  Since  t (p)  = - T/2irp  , we  have 

P V 

2 2 

Tp^(a)/Tp^(b)  =b  /a  in  any  equilibrium  configuration  of  the  body  ir- 
respective of  the  magnitude  of  the  applied  twist.  Thus  (3.23)  can  be 
written  as 


(a) 

(b) 


< 


T 

max 

T 

min 


The  torque -twist  relations  (3.11),  (3.16)  and  (3.  19)  are 
sketched  in  Fig.  2.  Clearly,  these  curves  are  anti- symmetric  with 
respect  to  the  axes. 

We  observe  from  the  preceding  calculations,  and  also  from 
Fig.  2,  that  we  have  not  as  yet  found  any  solutions  to  (2.  25)  - (2.  27) 
if  the  prescribed  twist  lies  in  one  of  the  intervals 


b 

J?F1 


H 


5F2 


> (3.24) 


J 


We  show  in  the  next  section  that  there  are,  in  fact,  no  smooth  solutions 
when  the  prescribed  twist  lies  in  these  ranges. 
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3 . 2 Non-existence  of  Smooth  Solutions 

We  now  establish  a sequence  of  lemmas  leading  to  a result 

which  is  in  fact  stronger  than  the  one  claimed  at  the  end  of  the  last 

section.  We  will  show  that  there  is  no  solution  0(r)  to  (2.25)  - (2.27) 

which  is  continuously  differentiable,  if  the  prescribed  twist  is  in  one 

of  the  intervals  defined  by  (3.24). 

Lemma  1;  There  is  no  continuously  differentiable 
solution  0(r)  to  the  differential  equation  (3.  1), 
where  T is  a constant  and  f is  a continuously  dif- 
ferentiable function  conforming  to  (2.28)  - (2.30), 
for  which  (3.2)  holds  at  some  radius  s , a<s<b. 

Proof;  Assume  that  there  exists  such  a solution  $(r)  and  suppose  that 

k(s)  = s0/(s)  = + k1  . (3.25) 

By  hypothesis  k(r)  = r$'(r)  is  continuous  on  (a,  b)  so  that,  in  particular, 
it  is  continuous  at  r = s.  Therefore,  given  any  number  s>0  , there 
exists  a number  6 = 6(*)>0  such  that  | k(s)  - k(r)  | < $ for  all  r such 
that  |r-s|<6(«).  Using  (3.25)  we  may  write  this  as 

|kj-k(r)|<c  for  all  |r-s|<6(e)  . (3.26) 

Recall  that  f(k)  has  a local  maximum  at  k = k^  , so  that  there 
is  a number  r)>0  such  that 

f(kj)if(k)  for  |kj-  k | <T|  . (3.27) 

Combining  (3.26)  with  (3.27),  we  have 

f(kj)if(k(r))  for  | r - s | < 6(T))  , 


which  on  using  (2.t>),  (3.  1)  and  (3. 25)  leads  to 


Note  from  (3.1)  and  (3.25)  that  T = -2trs  f(k.)  , whence  T<0 


Equation  (3.28)  now  requires  that 


which  is  impossible.  Consequently  there  cannot  exist  a solution  0(r) 
with  the  properties  we  assumed. 

The  cases  s0'(s)=-k.  , ± k?  can  be  dealt  with  similarly. 


Lemma  2;  Suppose  that  there  exists  a continuously 
differentiable  solution  0(r)  to  (2.25)  - (2.27), 
where  T is  a constant  and  f is  as  in  Lemma  1. 
Then 

(i)  |r$'(r)|<kj  on  (a,  b)  if  and  only  if  0^  is  in 
the  interval  (3.5). 

(ii)  kj  < | r0'(r)  | < k^  on  (a,  b)  if  and  only  if  0q  is 
in  the  interval  (3.29) 

(iii)  |r0'(r)|>k2  on  b)  if  and  only  if  0Q  is  in 
the  interval  (3.  18) 


Proof:  Considering  part  (i),  suppose  that  0^  is  in  the  interval  (3.  5) 
Byvirtueof  Lemma  1,  the  steps  leading  from  (3.5)  to  (3.8),  go 
through  even  when  0 is  merely  continuously  differentiable.  Thus 
necessarily  |r0'(r)|<kj  on  (a,b). 

Conversely,  suppose  that  |r0/(r)|<k.  on  (a,  b).  It  follows 


from  (2.25)  and  (3.25)  that 
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Proof:  This  result  follows  immediately  from  Lemmas  1 and  2.  For, 
suppose  that  there  is  such  a solution  0(r).  It  follows  from  Lemma  2 
that  we  must  have 

r0'(r)=±kj  or  ±1^  at  some  r,  a<r<b  . (3.36) 

But  Lemma  1 says  that  this  is  impossible. 

We  have  thus  shown  that  for  certain  ranges  of  the  prescribed 
twist,  there  is  no  solution  in  the  classical  sense  to  the  problem  under 


consideration. 
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4.  1 Weak  Formulation  of  Problem 

There  are  some  problems  of  considerable  physical  interest  in 
which  the  field  quantities  do  not  vary  smoothly  through  the  body.  Rice 


gives  some  examples  of  such  problems  in  [ 14].  We  have  observed  that 
the  problem  under  consideration  here  has  no  smooth  solution  for  certain 


ranges  of  the  applied  loading.  One  possibility,  which  we  shall  not  con- 
sider, is  that  the  tube  buckles,  possibly  into  some  unsymmetric  state 
of  plane  strain  at  such  a loading.  An  alternative  possibility  is  that  the 
tube  remains  in  a configuration  of  axisymmetric  plane  strain,  but  that 
now  the  field  quantities  are  no  longer  smooth  and  exhibit  certain  dis- 
continuities. This  latter  possibility  is  suggested  by  the  observation  in 
Section  2.2  that  the  displacement  equations  of  equilibrium  may  suffer 
a loss  of  ellipticity  at  certain  deformations  for  the  material  at  hand. 

In  particular,  in  view  of  known  results  in  the  theory  of  transonic  gas 
flows,  one  would  anticipate  the  occurrence  of  curves  across  which  the 
first  derivatives  of  the  displacement  field  suffer  jump  discontinuities, 
while  the  displacement  field  itself  remains  continuous. 

General  questions  concerning  the  possibility  of  the  change  of 
type  of  the  differential  equations  governing  finite  plane  elastostatics  and 
the  related  issue  of  the  existence  of  equilibirum  fields  with  discontinuous 
deformation  gradients  have  been  investigated  in  [4]  - [8].  Elastostatic 
fields  possessing  discontinuities  of  this  type  are  referred  to  as  "elasto- 
static shocks". 

We  now  relax  the  smoothness  demanded  of  the  local  twist  0(r) 
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and  the  pressure  field  q(r),  in  the  hope  that  this  will  enable  us  to  ex- 
plain what  happens  when  the  prescribed  twist  0^  is  in  one  of  the  inter- 
vals (3. 24). 

To  this  end,  let  7 be  a number  in  the  interval  [a,  b].  If  in  fact 
a<7<b  , we  will  now  require  that  0(r)  be  merely  twice  continuously 
differentiable  on  the  intervals  (a,  r)  and  (7,  b)  and  continuous  on  [a,  b]. 
The  stress  field  and  pressure  field  induced  by  the  deformation  (2.  1)  are 
to  be  continuously  differentiable  on  (a,7)  and  (7,  b)  while  the  traction  is 
presumed  to  be  continuous  at  r =7.  Accordingly,  we  have  admitted  the 
possibility  of  the  existence  of  a cylindrical  elastostatic  shock1  of  radius 
7 co-axial  with  the  cylindrical  region  ft  . 

The  global  balance  laws,  which  continue  to  be  meaningful,  now 
reduce  to  the  same  differential  equations  obtained  in  Section  2.  1 on 
(a, 7)  and  (7,b),  together  with  jump  conditions  at  r=7.  Accordingly 
we  now  have 


|?{r2f(r0/(r))}  = O , 


Z~W'(Z  + r20/2(r))  - 2r0'2(r)W/(2  + r20'2(r))  = |S  , 


(4.  1) 

on  a<r<7,  7<r<b, 

(4.2) 


instead  of  (2.20)  and  (2.21).  On  integrating  (4.  1)  we  have 


f(r0'(r))  = j on  (a,  r)  , f(r0'(r))  = - — on  (7,  b)  , (4.3) 

2*r  2ur 


where  T and  T are  (not  necessarily  equal)  constants.  Because  of 


We  formulate  the  problem  in  the  case  when  a single  elastostatic  shock 
exists.  We  will  find  that  this  suffices  for  our  purposes,  and  more  im- 
portantly, that  a configuration  involving  more  than  one  shock  is  neces- 
sarily unstable  (in  a sense  to  be  made  precise). 


4 
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(2.10),  (2.13),  (2.24)  and  (4.3)  we  have  that 

Tpt(r)  = '7X2  on  (»'*),  %*<*)  = on(r,b)  . 

2irr  Pf  2irr 

At  r = r,  equilibrium  considerations  require  that  the  tractions  be 
continuous.  Clearly,  this  is  equivalent  to 


(4.4) 


(4.5) 


TppF*,  = Tpp(Ff)  ' 


(4.6) 


Equations  (4.4)  and  (4.5)  lead  to 


- + 

T = T = T 


(4.7) 


We  therefore  have  the  following  problem  governing  the  local 
twist  0(r)  . Given  a number  0Q  , find  a function  0(r)  which  is  contin- 
uous on  [a,  b]  and  twice  continuously  differentiable  on  (a,r)  and 
(r,  b),  and  numbers  T,  r with  aSrsb,  such  that 


f(r0/(r))  = - 


2 on  a<r< r , r<r<b  , 


(4.8) 


0(a)  = 0Q  , 


(4.9) 


0(b)  = 0 


(4.  10) 


Integration  of  (4.2)  leads  to 
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m 


q(r)  = 


2W/(2  + rV2(r))  + j250/2(?)W/(2  + ?20/2(5))d5+qo  on  (a,r) 


2W/(2  + rV2(r))- j250'2(5)W/(2  + 520/2(5))d5+qo  on  (r,  b) 


(4. 11) 


which,  together  with  (2.  10)  and  (2.  13)  , gives 


t (r)  = 
PPV  ’ 


- J 250/2(5)W'(2  + 5 V2( 5))d5  - q. 


on  (a,  r)  , 


(4.12) 


J 2?0/2(5)W/(2  + 520/2(5))dS  - q0  on  (F,  b)  . 


(4.  13) 


We  see  from  (4.  12)  and  (4.  13)  that  the  traction  continuity  condition 
(4.6)  holds  if  and  only  if 


qo"qo  * 


(4.  14) 


Once  (4.8)  - (4.  10)  has  been  solved  for  the  function  0(r)  , equation 
(4.11)  together  with  (4.14)  directly  gives  the  pressure  field  q(r). 

4.  2 Weak  Solutions 

We  first  observe  that  the  Lemma  at  the  beginning  of  Section  3.  1 
continues  to  hold  if  we  replace  (a,b)  by  (rJf  r^)  where  rj  and  r2 
are  any  two  numbers  such  that  air^r^b  . This  result,  with  the 
particular  choices  r ^ = a , r^  = 7 and  r ^ = 7 , r = b , leads  to  the  con- 
clusion that  all  admissible  solutions  of  (4.8)  are  necessarily  such  that 
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r0'(r)  7*  ± kj  , ± k2  on  a<r<r,  r<r<b. 


(4.  15) 


Since  0'(r)  is  continuous  on  (a,  r)  it  now  follows  that  any  admissible 
solution  to  (4.8)  must  be  such  that  r0'(r)  takes  on  values  exclusively 
in  one  of  the  intervals  (-oo, -k2)  , (-k^-kj),  (-kj.kj),  (k^k-)  or 
(kj,  oo)  , at  all  points  in  (a, 7)  . The  same  must  be  true  on  (r,  b). 
Therefore  we  see,  because  of  (2.33)  and  (2.39)  , that  (4.8)  holds 


if  and  only  if 


r0'(r)  = 


-Fi(-L^)  on  (a,  r)  , 
2irr 


F.(— ^-y)  on  (r,  b)  , 

J Vtt  ^ 


(4.16) 


for  some  fixed  i,  j = 1,2,3  . 

Integrating  (4.16)  and  using  the  boundary  conditions  (4.9)  and 
(4.  10)  leads  to 


0(r)  = 


vJ?Fi(-7z)d!  °n  [a-r)  • 


“ ^.*>1  • 


(4.17) 


Finally,  we  require  that 


(4. 18) 
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since  the  local  twist  0(r)  given  by  (4.  17)  is  supposed  to  be  continuous 
at  r = r . 

Collecting  the  preceeding  results,  we  come  to  the  following  con- 
clusion. Given  a real  number  , if  there  exist  real  numbers  T and 
F , a^rsb  , such  that  (4.  18)  holds  for  some  fixed  choice  of  the  sub- 
scripts i,  j = 1,  2,  3 , then  (4.  17)  with  this  choice  of  T , F , i and  j 
is  a solution  to  (4.8)  - (4.  10)  at  the  given  . 

Clearly  in  the  case  when  i = j=  1,2,3  , (4.17)  and  (4.18)  de- 

scribe the  smooth  solutions  we  obtained  in  Section  3.1.  This  is  not 
surprising,  since  any  smooth  solution  of  (2.25)  - (2.27)  is  also  a 
solution  of  the  problem  in  its  weak  formulation.  Likewise,  in  the  parti- 
cular cases  when  F = a and  F = b,  (4.17)  and  (4.18)  are  readily  seen 
to  reduce  again  to  these  same  smooth  solutions.  A solution  defined  by 
(4.17)  and  (4.18)  is  therefore  not  smooth  only  if  i^j  and  a<F<b  . 

The  existence  of  a solution  (4.  17)  corresponding  to  the  pre- 
scribed value  of  the  twist  0q  is  contingent  upon  the  existence  of  numbers 
T and  F,  a^F^b  , such  that  (4.18)  holds.  We  now  examine  this  lat- 
ter issue.  First  note  that  since  (4.  18)  furnishes  only  one  scalar 
restriction  on  the  two  numbers  T and  r"  , we  expect  that  if  there  are 
values  of  T and  F conforming  to  (4.  18),  then  there  would  in  fact 
be  many  such  values.  If,  therefore,  we  momentarily  imagine  specify- 
ing both  and  T , we  may  pose  the  following  question:  at  each  fixed 
choice  of  the  subscripts  i,  j=  1,2,3  , ij^j  , for  what  values  of  the  pair 
(0q,  T)  will  (4.18)  determine  a value  for  r , a^F^b?  We  will,  with 
no  loss  of  generality,  restrict  attention  to  the  first  quadrant  of  the 
(0q,  T)-plane.  We  now  show  that  for  each  fixed  choice  of  the  subscripts 
i,j  = 1,2,3  , i^j  , there  is  a simply  connected  closed  region  in  the 


-32- 


first  quadrant  of  the  (0Q-  T)  plane  such  that  (4.  18)  determines  a value 
for  7 if  and  only  if  (0q,T)  is  in  A..  . Furthermore,  this  value  of  r 


is  unique. 


To  this  end,  define  the  functions  0^  , i,j  = l,2,  3 , i^j  , by 


^•T>  = J?Fi(^)d5+J?Fj(^p>5  Bij  • W.19) 

^ r 

where  the  domains  of  definition  B. . of  the  functions  0..  on  the  (r,  T)- 

ij  ij  ' 7 

plane  are  given  by 


3 1= 

{(r.T)| 

as 

7s 

b , 

0 -2 

2irr  T 

13= 

{CF,  T)| 

as 

rs 

b , 

2irb^T 

l 

2 1= 

{(r,  T)  | 

as 

7^ 

b , 

2it72t] 

12= 

{C?,T)| 

a s 

7s 

b , 

2iTb2T 

l32  = 

{ Cr , T)  | 

as 

r s 

b , 

2trb2T 

23= 

{(r,X)| 

as 

7s 

b , 

2TTb2T^ 

(4.  20)J 


We  now  consider  the  case  i = 3 , j=  1 in  detail.  For  each  fixed 

value  of  7 in  [a,b]  , it  follows  from  (4.19)  , (4.20)  that  0q=0^j(t,  T) 

defines  a segment  of  a smooth  curve  on  the  (0q,  T)-plane  for 
—2  —2 

2irr  t . sT£2irr  t . Therefore,  we  have  a family  of  such  curves 
min  max  ’ 7 

on  the  (0q  ,T)-plane,  each  corresponding  to  a different  value  of  r in 


Because  of  (2.32)  , (2.33)  one  sees  that  these  are  the  largest  possible 
domains  of  definition  of  the  functions  0^j  . In  the  case  when 
a^Tmax<*)^Tmin  ~ 80  (3-23)  does  riot  hold  — we  see  that  B13  , B12 

and  B23  are  empty.  In  this  case,  therefore,  solutions  (4.  17)  with 
(i,  j)  = (1,  3)  , (1,  2)  and  (2,  3)  do  not  exist. 


[a,  b]  , and  all  of  them  having  their  end  points  on  the  curves 


^0=^3i^T^2irTmin'  T>  and  ^0=  ^3  l^T/^2lTTmax'  * Since  by  <2*38)  411(5 
(4. 19)  we  have 


803l(r.  T) 
07 


on  B31  , 


(4.21) 


it  follows  that  the  different  members  of  this  family  of  curves  do  not  in- 
tersect each  other.  Furthermore,  a curve  corresponding  to  a larger 
value  of  7 lies  to  the  right  of  a curve  corresponding  to  a smaller  value 
of  r . And  finally,  since  03j  depends  continuously  on  7 , these  curves 
span  a simply  connected  region,  , in  the  (0q  , T)-plane.  From  the 

above  discussion  it  follows  that  A3j  is  the  closed  region  bounded  by 
the  curves  0Q=031(a,T),  0Q=031(b,T),  0Q=  03 i(VT/2TrTmin»  T)  411(1 
V*31<^W.  T)  ’ A sketch  of  this  regie  , together  with  the  span- 
ning family  of  curves,  is  shown  in  Fig.3(i)  . The  fact  that  a curve  cor- 
responding to  a larger  value  of  7 is  to  the  right  of  a curve  associated 
with  a smaller  value  of  7 is  indicated  in  Fig.3(i)  by  the  arrow  labelled 
"direction  of  increasing  7".  Since  there  is  exactly  one  of  these  curves 
passing  through  any  point  in  A3j  , it  follows  that  there  is  a unique 
number  r associated  with  every  point  (0q,  T)  in  A3j  , such  that 
0O=03i(r,  T).  This  is  what  we  set  out  to  establish.  We  may  express 
this  analytically  as  follows:  there  exists  a function  ?3j  , defined  on  A3j, 
such  that  r determined  by 

7 = ?31(0o,T)  (4.22) 


31  ‘ 


conforms  to  0q  = 03j(7,  T),  i.e.  0Q=  ^3 1(^3 1^0’  T^’  T)  on  A 
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Summarizing  the  results  for  this  case,  we  have  that,  if  and 
T are  numbers  such  that  (0A,  T)  is  in  A (the  region  PQRS  in 

u 31  ' 

Fig.3(i))#  then  there  is  a unique  number  r (a^rsb)  such  that  (4.  18) 
holds  (with  i = 3,  j = 1).  Equation  (4.  17)  , with  these  values  of  T,  r,  i 
and  j , is  a solution  to  (4.8)  - (4.  10)  at  that  value  of  . 

The  other  cases  — corresponding  to  the  remaining  choices  of  the 
subscripts  i,  j — may  be  likewise  examined.  In  each  case  we  find  a 
simply  connected  closed  region  , shown  in  Figs.  3-5,  such  that,  if 
and  T are  numbers  with  (0q,  T)  in  A^  , then  there  is  a unique 
number  "r  (a^7^b)  such  that  (4.  18)  holds  for  that  choice  of  i,j  . 
Equation  (4.  17)  then  provides  the  corresponding  solution  0(r).  Accord- 
ingly, in  each  case  there  exist  functions  r^.  defined  on  A^  , such  that 


r = ?..(VT>  on  A.. 


(4.23) 


conforms  to  0Q=  0^(7,  T)  . 

The  composite  torque-twist  diagram,  wherein  all  of  these  admis- 
sible regions  A.j  together  with  the  torque-twist  curves  for  the  smooth 
solutions  are  sketched  on  one  figure, is  shown  in  Fig. 6.  We  observe  that 
the  admissible  regions  A.,  "fit"  appropriately  between  the  torque-twist 
curves  associated  with  the  smooth  solutions  (Fig. 2).  Therefore  cor- 
responding to  any  given  value  of  the  twist  0q  we  now  have  a solution. 
However,  we  are  now  faced  with  the  unsatisfactory  situation  in  which 
there  is  an  infinite  number  of  admissible  solutions  at  certain  values  of 
the  prescribed  twist  0q  . 

We  observe  from  Fig. 6 that  at  sufficiently  small  twists  0Q 
(s  0 ) and  at  sufficiently  large  twists  0q  (a  0q)  we  have  a unique 
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solution,  which  is  smooth.  When  the  prescribed  twist  0^  is  in  one  of 
the  intervals  * ^N<^0<^R  * we  ^ave  411  number  of 

solutions,  all  of  which  are  weak  solutions.  In  the  remaining  intervals, 
we  have  an  infinite  number  of  solutions,  one  of  which  is  smooth,  all  the 
rest  being  weak  solutions. 

Observe  from  Figure  6 that  even  a knowledge  of  both  0q  and  T 
may  be  insufficient,  in  some  cases,  to  determine  a unique  solution.  For 
example,  there  are  four  solutions  corresponding  to  any  point  in  PMNK, 
one  for  each  of  the  pairs,  (i,  j).=  (1,  2)  , (2,1),  (1,3)  and  (3,  1).  We 
remark  that  at  any  point  (0q,  T)  on  PS  or  RQ  there  is  in  fact  only 
one  solution  — the  smooth  one.  One  sees  this  from  (4.  17)  , (4.  18)  , 
since  all  of  the  weak  solutions  at  such  a point  have  either  7 = a or  7 = b 
(see  Figs.  3-5).  Likewise,  at  any  point  on  MN  we  only  have  smooth 
Solution  2 or  the  weak  solutions  (i,  j)  = (1,3)  , (3,  1). 

Finally,  we  observe  that  it  is  convenient  to  visualize  the  various 
solutions  as  follows.  Consider  for  example  a weak  solution  with  i = 3 , 
j=  1.  Let  7 denote  the  radius  of  the  associated  shock.  Let  A,  B,  C 
and  D be  points  on  a radial  line  in  the  cross-section  II  of  the  tube  in 
the  undeformed  configuration,  see  Fig.  7(a),  such  that  A and  D are 
at  the  inside  and  outside  boundaries  respectively,  while  B and  C are 
points  just  inside  and  outside  the  shock-line.  The  solution  at  hand  is 
given  by  (4.17)  with  i = 3 , j = 1 . If  we  use  this  to  compute  r0'(r)  and 
then  plot  the  points  with  coordinates  (|r0'(r)|  , T/2rrr  ) (suppose  T>0) 
for  each  r in  the  intervals  »sr<7  , T<rSb  , we  obtain  the  curves 
AjBj  and  CjDj  (typically)  shown  in  Fig.  7(b).  The  graph  of  f(k)  has 
been  superimposed  on  this  diagram.  The  abscissa  of  any  point  on  AjBj 
or  CjDj  gives  the  value  of  the  local  amount  of  shear  |r0'(r)|  at  the 


m 
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corresponding  point  in  the  tube,  while  the  ordinate  gives  the  correspond- 
ing shear  stress  |t^  | . Observe  from  Fig.  7(b)  how  the  local  amount 
of  shear  varies  continuously  on  either  side  of  the  shock  but  suffers  a 
jump  discontinuity  across  it.  The  shear  stress  , on  the  other  hand, 

is  seen  to  vary  continuously  throughout  the  tube.  If  we  refer  to  the  por- 
tions of  the  curve  f(k)  vs.k  between  Oskskj  , k^sksk^  and 
k2Sk<co  as  the  first,  second  and  third  branches  of  f(k)  respectively, 
we  see  that  this  solution  (i  = 3,  j=  1)  is  associated  with  the  third  and 
first  branches  of  f(k)  , with  the  region  inside  the  shock-line  associated 
with  the  former  branch.  In  general,  the  weak  solution  (i.j)  is  associ- 
ated with  the  i*^  and  j**1  branches  of  f(k).  with  the  part  of  the  tube 
inside  the  shock-line  corresponding  to  the  i^1  branch. 

We  see  from  this  and  Section  2.  2 that  the  type  of  these  weak 
solutions  is  mixed,  in  general.  The  displacement  equations  of  equili- 
brium are  elliptic  on  that  part  of  n for  which  a<r<T  and  non-elliptic 
where  r<r<b  , at  solutions  with  (i.j)  = (1,2)  , (3,2)  , while  they  are 
elliptic  where  r<r<b  and  non-elliptic  where  a<r<r,  at  the  solutions 
(i*  j)=  (2,  1)  , (2,3).  In  the  case  of  the  solutions  corresponding  to 
(i*j)  = (1*3)  * (3,  1)  , the  displacement  equations  of  equilibrium  are  el- 
liptic everywhere  in  II  where  r^T. 


1 


5.  1 Dissipativity 

The  lack  of  uniqueness  encountered  in  the  preceeding  section  is 
not  unexpected  since  we  had  enlarged  the  admissible  class  of  solutions 
there.  In  such  circumstances,  it  is  usually  the  case  that  not  all  of  the 
solutions  admitted  by  the  differential  equations  are  physically  reasonable. 
In  gas  dynamics,  for  example,  there  are  problems  in  which  the  differen- 
tial equations  admit  solutions  which  are  unacceptable  since  they  are  as- 
sociated with  a decrease  of  entropy.  It  is  essential  therefore  to  intro- 
duce additional  criteria  which  will  single  out  a physically  admissible 
solution. 

Knowles  and  Sternberg  proposed  such  a criterion  in  [ 6 ] , in  the 
context  of  finite  elastostatics,  which  they  referred  to  as  the  dissipativity 
inequality.  A thermodynamic  motivation  for  this  inequality,  stemming 
from  the  Clausius- Duhem  inequality,  was  given  by  Knowles  in  [ 7 ] . The 
dissipativity  inequality  is  essentially  an  expression  of  the  physically 
reasonable  idea  that  the  rate  at  which  elastic  energy  is  beina  stored  in 
any  part  of  the  body,  in  some  quasi-static  process,  cannot  exceed  the 
rate  at  which  work  is  being  done  on  that  part. 

We  now  examine  the  implications  of  the  dissipativity  inequality 
in  the  context  of  the  present  problem.  While  we  could  specialize  the 
general  dissipativity  inequality  given  in  [7]  for  our  problem,  it  is  il- 
. lustrative  (and  equally  easy)  to  derive  it  from  first  principles. 

We  now  consider  a quasi-static  time-dependent  family  of  equili- 


brium solutions.  The  time  t merely  plays  the  role  of  a history 
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parameter  and  no  inertia  effects  are  accounted  for.  Accordingly,  we 
are  concerned  with  a one-parameter  family  of  functions  0(r,  t)  , depend- 
ing on  the  parameter  t in  some  time  interval  T , such  that  at  each  t 
in  T , 0(r,t)  is  a solution  to  (4.8)  - (4.  10)  . The  torque,  twist  and 

shock  radius  are  all  time  dependent  now,  and  we  write  T(t)  , 0g(t)  and 
'r(t).  It  is  convenient  to  set 


II(t)  = { r |‘r(t)<  r<  b}  , 
fi(t)  = [r  |a<r<r(t)}  , 


f o r t in  T . 


(5.1) 


Then  0(»,t)  is  continuous  on  [a,  b]  and  twice  continuously  differentiable 

+ 

on  II  and  n at  each  t in  I . Furthermore 


f(r0r(r,t)) 


= - ^^2  on  ^(t)  and  fl(t) 


(5.2) 


0(a,t)  = 0o(t)  , 


(5.3) 


0(b,t)  = O , 


(5.4) 


at  each  t in  I , Here  0q (t ) is  the  prescribed  twist,  and  we  suppose 
it  to  be  continuous  and  piecewise  continuously  differentiable  on  T . In 
certain  discussions,  as  we  observed  previously,  it  will  be  temporarily 
necessary  to  imagine  that  T(t)  is  also  specified.  In  such  circumstances, 
we  presume  T(t)  to  possess  the  same  smoothness  as  0p(t)  on  T. 

It  is  convenient  to  set 


k(r,t)=  r0r(r,t)  , 


*We  use  the  notation  0 _9£(r,t)  0 t 

r or  t St 
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k(t)  = r(t)0r(r(t)+,t)  , 


(5.6) 


k(t)  = r(t)0r(r(t)-,t)  , (5.7) 

with  the  understanding  that  when  T(t)  = b we  take1  r(t)+  = b and  when 
TCt)  = a we  take  T(t)-  = a . k and  k represent  the  instantaneous  local 
amounts  of  shear  at  points  just  outside  and  inside  the  shock,  respectively. 

We  will  now  require  that  at  each  instant  in  I , the  rate  at  which 
the  external  forces  on  the  tube  are  doing  work  should  not  be  less  than 
the  rate  of  increase  of  the  stored  energy,  i.e.  we  demand  that 


W(2  + k (r,  t))2irrdr  for  all  t in  T 


(5.8) 


We  may  evaluate  the  right  hand  side  of  (5.  8),  using  (5.  6)  and  (5.7), 
as  follows: 


b 

^ J W(2  + k2(r,  t))2irrdr 
a 

r(t)  b 

= |^Jw(2  + k2(r,t))2irrdr+^  J W(2  + k2(r,  t))2irrdr  , 
a r(t) 

b 

= {W(2  + k2)-W(2  + £2)}  2ifr(t)~(t)  + J 4rrrk  Jj^  W#(2  + k2)dr  . (5.9) 


Using  (2.24),  (5.2)  and  (5.5)  in  (5.9),  gives 


^his  is  admissible  since  we  observed  in  Section  4.  2 that  when  r = a 
or  b , the  solution  is  in  fact  smooth.  Thus  exists  there. 
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^Jw(2  + k2)2Trrdr 


= { W(2  + k2)  - W(2  + k2)32nr(t)|^  - T(t)  J 0rt(r,  t)dr 


However,  because  of  (5.3)  and  (5.4),  we  have 


a# 

J *rt(r.  t)dr  = 0t(b,  t)  - 0t(r(t)+,  t)  + 0t(r(t)-,  t)  - 0t(a,  t)  , 


= ‘dt0O(t)"#t(F(t)+'t)  + 0tfF(t)"»t)  ' 


so  that  we  may  write  (5.  10)  as 


^ J W(2  + k2)2irrdr  = f W(2  + k2)  - W(2  + it2)}  2nr(t)^f+  T(t )|^  0Q(t) 


+ T(t){0t(F(t)+,t)-0t(7(t)-,t)} 


Since  the  displacements  are  continuous  across  the  shock,  we  have 


0(r(t)+,t)  = 0(r(t)-,t)  for  t in  T , 


which  when  differentiated  with  respect  to  t leads  to 


0rCF(t)+.  t)ff  (t)  + 0t(r(t)+,  t)  = 0r(r(t)-,  t)|f(t)  + 0t(r(t)-,  t) 


(5.  10) 


(5.11) 


(5.12) 


(5.13) 


(5.14) 


Using  (5.  6),  (5.  7)  and  (5.  14)  in  (5.  12)  gives 
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|t  J W(2  + k2)2irrdr  = { W(2  + k2)  - W(2  + k)2}2Trr|f  + T(t)|^  0Q(t) 
a 

+=ff)T(t){k(t)-k(t)}|f(t)  , (5.15) 


which  because  of  (5.2),  (5.6)  can  be  written  as 


J W(2  + k2)2irrdr  = T(t)^0Q(t)  + { W(2  + k2)  - W(2  + k2)}21r7|f 


+ 27t7f(k)(k-k)^(t)  . 


(5.16) 


The  dissipativity  requirement  (5.8)  can  now  be  written  as 


{ W(2  + k2)  - W(2  + k2)  - f (k)(k  - k)}  2irT~i  0 on  T , 


or  alternatively,  because  of  (2.24),  as 
+ 


| Jf(5)d5-(k-k 


k)f(it) 


2irr^^0  for  all  t in  T . 


(5.17) 


This  is  the  form  of  the  dissipativity  inequality  that  we  shall  find  useful. 

It  follows  from  the  results  of  Section  4.2  that  all  admissible 
quasi-static  families  of  equilibrium  solutions  are  of  the  form 


1 Note  from  (5.2)  that  f(k)=f(ic)  . 
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r 


0O(t)-J?Fi(7^)«  0“  fi(t)  , 

a ^tt5 


0(r#  t)  =< 


J,?Fj(f^)dS  on  ^(t)  ' 

r 


subject  to  the  restriction 


(5.  18) 


r(t)  b 


for  i,  j = 1,2,3  and  for  t in  I , We  now  proceed  to  apply  the 
dissipativity  inequality  (5.  17)  to  the  various  families  of  solutions  rep- 
resented by  (5.18),  (5.19). 

We  first  note  that,  if  at  some  instant  t we  have  a smooth  solu- 
tion, then  (5. 17)  holds  at  that  instant  by  virtue  of  the  continuity  of  0 
i.e.  since  k = £ . Therefore,  we  may  restrict  attention  to  the  cases 
for  which  i^  j in  (5. 18)  , (5.  19)  , and  to  times  in  I for  which 


a<r(t)<b  . 


Equations  (5.6),  (5.7)  and  (5.18)  now  give 


k(t)  = -F 


.(-If' 

1 srrrft) 


See  discussion  following  (4.  18). 


(5.20)1 


(5.21) 


rry 
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so  that  we  may  use  (2.28),  (2.39),  (2.40),  (2.41),  (5.2),  (5.6)  and 
(5.  21)  to  write  (5.  17)  as 


F.h(t)] 

2irF(t)  £ (t)J  ff(5)  - n(t)  } dS  * 0 , Ti(t)  = 

F.[t)(t)] 


1M 

2ir7^(t) 


(5.22) 


for  all  t in  T for  which  (5.20)  holds. 

It  is  convenient  to  define  the  functions  Aj  and  by 

F2<t) 


Ai(t)=J  {f(5)  - t}d5  for 


max 


(5.23) 


Fj(t) 


F3(t) 


A2(t)=J  {t -«{)}«  for 


max 


(5.24) 


F2(t) 


These  functions  Aj  and  A2  have  the  following  geometrical  interpreta- 
tion. If  in  Fig.  1 we  draw  a line  parallel  to  the  k-axis  at  a distance  t 

above  it  (t  . st^t  ) , then  A.(t)  and  A,(t)  are  the  areas  of  the 
min  max  i <■ 

two  loops  formed.  Clearly 


A.(t  ) = A_(t  , ) = 0 , 
lv  max'  2 min 


(5.25) 


A.(t)>0  for  t . £t<t 

1'  ' min  max 


(5.26) 


A,(t)>0  for  t . <t£t 
2'  ' min  max 


(5.27) 
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Using  (5.23)  and  (5.24)  in  (5.22)  leads  to 


-.»)•*, W) 

2irr  2irr 


2wrd?i0  for  <U>=(3.1)  , 


2vr^0  for  (i,j)  = (l,3)  , 


■Al(“^i)2irT|fa0  for  (iJ)=(2,  1)  . 

2ttt 

Ai(~^ii)2ir'?lf s 0 for  (i*j)  = (1*2)  » 

2irr 


A2(“^li)2lTrft'S:0  for  J)  = <3*  2>  * 

2irr 


A2(“^Z2)2lT7|fi0  for  (i»j)  = (2.3)  . 

2trr 


(5.28) 

(5.29) 

(5.30) 

(5.31) 

(5.32) 

(5.33) 


in  each  of  the  different  cases. 

Now  consider,  for  example,  the  case  (i,j)=(2,l),  i.e.  suppose 
that  for  all  times  sufficiently  close  to  some  tj  in  I,  the  quasi-static 
family  of  solutions  (5.18)  has  i = 2,  j=  1.  We  then  have  from  (2.32), 

(2.33),  (5.  1),  (5.  18),  (5.26)  and  (5.30)  that  the  dissipativity  ine- 
quality is  satisfied  at  a time  tj  for  which  (5.20)  holds  if  and  only  if 

(5.34) 

I 1 

As  previously  observed,  in  the  event  that  (5.20)does  not  hold,  so  that 

(5.35) 


r(tj)  = a or  b , 
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the  dissipativity  inequality  holds  without  need  for  any  restrictions 
such  as  (5.34).  The  meaning  of  these  restrictions  is  most  transparent 
when  viewed  in  the  torque-twist  diagram  (see  Fig.  8).  With  no  loss  of 
generality  we  suppose  that  T(t)  and  0g(t)  are  non-negative  for  all  times 
in  T . We  shall  refer  to  the  piecewise  smooth  oriented*  curve  T in  the 
torque-twist  plane  defined  by  ^q=#q(*)  * T = T(t)  for  t in  T as  the 
loading  path.  By  hypothesis,  for  all  values  of  t sufficiently  close  to 
tj  , the  loading  path  T lies  in  . Let  Z = (0Q(tj)  , T(tj))  be  the 

point  on  T corresponding  to  t = tj  . 

2 

Recall  that  the  region  A^ j is  spanned  by  a one-parameter 
family  of  curves  0q=02j(t,  T)  , a^rsb  , and  that  a member  of  this 
family  of  curves  corresponding  to  a larger  value  of  the  parameter  "r 
lies  to  the  right  of  a curve  corresponding  to  a smaller  value.  Let  C be 
the  particular  member  of  this  family  with  equation  0^  = 0-,  j(r(t j)  , T)  , 
so  that  C passes  through  Z , (see  Fig.  8).  It  follows  that  the  shock 
radius  7 corresponding  to  any  point  in  A^j  to  the  right  of  C is  greater 
than  'r(tj)  , while  at  a point  to  the  left  of  C , it  is  less  than  "r(tj)  . 
Therefore,  dissipativity  — (5.34)  — requires  that  the  loading  path  T 
should  be  oriented  at  Z in  such  a way  that  it  does  not  point  to  the  right 
of  C , provided  Z is  not  a point  on  PS  or  MN.  This  is  shown  in 
Fig.4(i)  as  well,  in  which  the  cluster  of  arrows  indicates 
the  admissible  orientations  of  a loading  path  through  a typical  point. 

This  is  true  for  all  points  in  A£j  except  for  those  which  lie  on  PS 
and  MN.  At  a point  on  these  curves  the  loading  path  may  be  arbitrarily 

*r  is  oriented  in  the  direction  of  increasing  time. 

2See  Fig.4(i). 
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oriented,  by  virtue  of  (5. 35). 

Clearly  we  can  analyze  the  other  cases  in  an  entirely  analogous 
manner.  We  find  that  dissipativity  is  essentially  equivalent  to 


3H 

if 

(i.j)  = 

= d,2)  , 

(5.36) 

if 

U.j)  = 

= (3,2)  , 

(5.37) 

if 

(i,  j)  = 

= (2,3)  , 

(5.38) 

and  these  are  geometrically  interpreted  in  Figs.  4 and  5 as  before.  The 
only  exceptions  to  (5.36)  - (5.38)  are  respectively  at  points  on  the 
curves  PK  , MN  and  MN,  RQ  and  MN  , LR  , whereat  the  orientation 
is  arbitrary. 

Equations  (5.28)  and  (5.29)  — i.e.  the  cases  (i,  j)  = (3,  1)  and 
(1,3)  — can  also  be  similarly  examined,  taking  care  now  to  note  that 
[ A^(t)  - A j(t))  is  not  always  of  the  same  sign.  If  we  set 

A(t)  = A2(t)-A1(t)  for  TminSTSTmax  • <5'39> 


where  Aj  and  A 2 are  as  defined  previously,  we  find  because  of  (5.23) 
(5.27)  and  (5.39)  that 


A(t  )<0  , A(t  )>0  , 
' mm  ' ' max'  ’ 


3^(t)>0  for  t . stSt 
dT  min  max 


(5.40) 


(5.41) 


Since  A(t)  is  continuous,  it  follows  from  (5.40)  and  (5.41)  that  there 

is  a unique  number  T in  (t  , t ) such  that 

c min  max 


l 


— 
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A(tc)  = 0 , 


(5.42) 


A(t )>0  on  (Tc,Tmax]  , A(t)< 0 on  [Tmin,Tc)  . (5.43) 


The  number  r is  shown  in  Fig.  1;  since  A.(t  ) = A,(t  ) , the 

C * C c C 

two  hatched  regions  are  of  equal  area.  The  dissipativity  conditions 
(5.28)  and  (5.29),  because  of  (5.39)  , (5.42)  and  (5.43)  , are  equiva- 

lent to 


dt  2„r2  c 


i£i0  if-llU*  . 

dt  - —2  c ’ 

2nr 


— is  arbitrary  if  -1 — W = t 
dt  2nr2  c 


(i,j)=(3.1)  , 


(5.44) 


if4>r 
dt  2*?  C 


>T„  » 


^£*0  if  , 

dt  _ _2  c ' 

2irr 


is  arbitrary  if  = t , 


(i. j)  = a,  3)  . 


(5.45) 


Consider  the  case  (i,  j)=  (3,  1)  . One  shows  easily  that 
*O=031(,/T/2itVT)  iB  a curve  in  A^  ^ which  qualitatively  looks  as  shown 
in  Fig.3(i)—  curve  XY.  Corresponding  to  any  point  on  this  curve. 
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we  have  T/2tt72=tc  while  it  is  readily  seen  that  at  any  point  above  or  below  XY 
T/2irr2  is  respectively,  greater  than  or  less  than  tc  . Accordingly,  the 
dissipativity  inequality  is  equivalent  to  the  first,  second  and  third  of 
(5.44)  at  points  in  A3J  , respectively,  above,  below  and  on  the  curve 
XY.  The  arrows  in  Fig.3(i)  indicate  the  admissible  orientations  of  a 
loading  path  at  some  typical  points  in  Ajj  . As  before,  the  orientation 
at  points  on  PS  and  QR  is  arbitrary.  The  solution  (i,j)  = (l,3)  may 
be  similarly  interpreted,  as  shown  * in  Fig.3(ii). 

5.  2 Consequences  of  Dissipativity 

The  dissipativity  inequality  was  introduced  in  the  hope  that  it 
would  single  out  a physically  admissible  solution  from  among  the  many 
available  equilibrium  solutions.  We  now  demonstrate  that,  if  we  re- 
quire the  local  twist  0(r,  •)  to  be  continuous2  on  I at  each  r in 
[a»h]  # and  if  we  suppose  that  the  body  was  in  an  undeformed  config- 
uration at  some  time,  then  a configuration  corresponding  to  solutions 
(i»  j)  = (1>  2),  (2,  1),  (2,  3)  (3,  2)  or  smooth  Solution  2 cannot  be  attained 
at  any  subsequent  time. 

First,  omit  the  weak  solutions  (i,j)  = (1,  3)  and  (3,1)  from  dis- 
cussion. We  observe  from  Fig.  4 that  any  loading  path  in  Fig.  6 con- 
forming with  the  dissipativity  inequality  and  starting  from  O is  neces- 
sarily confined  to  the  curve  OP  for  all  subsequent  time.  Note  similarly, 

1An  examination  of  the  details  of  the  curve  0O=  0l3(n/t/2ttto  , T)  show 

that  it  is  possible  for  this  curve,  depending  on  the  specific  geometry 

and  constitutive  law,  to  intersect  a different  pair  of  boundaries  of  A13 

than  shown  in  Fig.  3(ii).  The  figure  is  drawn  for  ba/aaT  . <T<a  /baT 

min  c max 

2 

Note  that  despite  the  presumed  continuity  of  0<j(t)  and  T(t),  0(r,*) de- 
fined by  (5.  18)  is  not  necessarily  continuous  on  T , since  the  subscripts 
i and  j may  change  values  at  certain  times. 


A 


-49- 


from  Pigs.  5 and  6,  that  any  admissible  loading  path  starting  from  O7 
is  likewise  restricted  to  O'R  for  all  subsequent  time.  The  only  possi- 
ble way  of  achieving  a solution  (i,  j)  = (1,  2),  (2,3),  (3,2),  (2,1)  or 
Solution  2 is  then,  by  virtue  of  a loading  path  which  is  associated  with 
one  of  the  solutions  (i,  j)  = (3,  1),  (1,3)  for  some  time  interval  less  than 
some  time  tj  , and  with  one  of  these  solutions  after  time  tj  . One  sees 
readily  from  (5.18)  however,  that  this  involves  a discontinuity  in  0(r,  • ) 
at  the  time  tj  . Since  we  have  disallowed  this  possibility,  we  now  con- 
clude that  a configuration  corresponding  to  any  solution  associated  with 
the  second  branch  of  the  graph  of  f(k)  vs.  k cannot  be  attained  through 
a dissipative  quasi-static  deformation  process.  These  are,  incidentally, 
the  solutions  at  which  the  displacement  equations  of  equilibrium  are  non- 
elliptic  somewhere  in  n . 

However,  even  if  we  now  discard  the  solutions  associated  with 
the  second  branch  of  f , we  would  not  have  overcome  our  troubles  with 
non- uniqueness.  For  example,  consider  the  Solutions  1,  3 and 
(i»j)  = (3,  1).  The  appropriate  torque-twist  diagram  is  shown  in  Fig.  9. 

If  we  imagine  gradually  increasing  the  applied  twist  0Q  from  zero,  the 
only  available  loading  path  initially  is  OS.  During  the  next  stage, 

' dissiPativity  - see  Fig.3(i)-  disallows  all  loading  paths 
except  SX.  Once  the  applied  twist  0q  exceeds  the  value  0^  , however, 
the  loading  path  could  lie  anywhere  in  PQYX  , and  we  have  no  criterion 
for  deciding  which  path  to  follow.  Eventually,  for  0Q>0Q  » we  are 
restricted  to  the  path  QO'.  Likewise,  during  a steady  decrease 
of  the  applied  twist  the  loading  path  would  be  restricted  to  O'QY 
then  allowed  to  follow  an  arbitrary  path  (consistent  with  dissipativity) 
in  XYRS  and  finally  restricted  to  SO.  It  is  interesting  to  note 
that  if  in  either  case  the  loading  path  lies  on  the  curve  OXYO' 


i 
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then  the  quasi-static  process  is  dissipationless  in  the  sense  that  (5.8) 
would  hold  with  equality  at  every  instant  t . 

It  is  therefore  imperative  that  we  seek  an  additional  — or  possibly 
an  alternative  — physical  criterion,  to  the  dissipativity  inequality,  that 
would  sort  out  more  completely  the  issue  of  non -uniqueness. 


-51- 


6.  1 Preliminaries  on  Stability 

We  now  look  into  the  possibility  of  using  a stability  criterion 
instead  of  the  dissipativity  inequality,  in  order  to  single  out  a physically 
admissible  solution  to  the  boundary-value  problem  under  consideration. 

We  draw  attention  to  the  fact  that  we  will  not  make  use  of  the  partial 
success  achieved  through  the  dissipativity  inequality,  since  we  are 
at  present  examining  the  possibility  of  an  alternative  — rather  than 
additional  — criterion. 

The  notion  of  stability  that  we  will  use  is  a static  one  based  on 
the  energy  criterion.  * According  to  this,  an  equilibrium  configuration 
of  a body  is  stable  if  the  work  done  by  the  external  loads  in  every 
sufficiently  small  kinematically  possible  virtual  displacement  from  this 
equilibrium  configuration  is  less  than  the  corresponding  increase  in  the 
stored  energy. 

In  order  to  mathematically  formulate  this  criterion,  it  is  necessary 
to  decide  on  a measure  for  the  virtual  displacements  and  to  specify  the 
behavior  of  the  applied  loading  during  a virtual  displacement.  We  will 
consider  two  possibilities  — stability  under  dead  loading  and  stability 
with  fixed  boundaries. 

First  consider  the  case  in  which  we  have  dead  loading  on  the 
inner  surface  of  the  tube  while  the  outer  surface  remains  fixed.  The 
torque  T then  remains  constant  during  a virtual  displacement.  Let 


See  page  195  of  [ 1 5 ] for  a discussion  of  this  criterion. 
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0(r)  be  the  equilibrium  solution  whose  stability  we  wish  to  investigate. 
Define  the  potential  energy  functional  V { \|r } by 

b 

V{*}  = J(W(2  + rV2)  - W(2  + rV2)32Tirdr  - T{*(a)  - 0(a)}  (6.  1) 

a 

for  all  functions  i(r(r)  in  some  set  X • In  order  to  interpret  f(r)  as  a 
virtual  twist*  from  the  undeformed  configuration,  we  suppose  that  X 
is  the  set  of  all  functions  which  are  continuous  and  have  piecewise  con- 
tinuous first  and  second  derivatives  on  [a,  b],  and  for  which  t(b)  = 0. 
Since  this  limited  degree  of  smoothness  is  all  that  is  required  of  an 
equilibrium  solution  0(r)  , it  seems  reasonable  not  to  impose  more 
severe  smoothness  requirements  on  the  virtual  displacement.  Finally, 
we  will  measure  the  departure  of  a virtual  twist  from  the  solution 
0 by 

b 

lit  - 0 1 1 =j2Trr3[f/-0/]2dr  . (6.2) 

a 

We  now  say  that  an  equilibrium  solution  0(r)  is  s table  if  there 
exists  some  number  «>0  such  that 

V{f}>0  for  all  functions  t in  X (6.3) 

for  which  1 1 ♦ - 0 1 1 < « » t ^ 0 • 

A solution  which  is  not  stable  is  unstable,  i.e.  if,  given  any  number 
«>0  , there  exists  a function  ^(r)  such  that 


We  restrict  attention  to  purely  circumferential  virtual  displacements. 
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V{*}S0  , ||*  - 0||<C  . ♦ , (6.4) 

the  equilibrium  solution  0(r)  is  unstable. 

We  now  determine  from  (6.  1),  (6.3)  a sufficient  condition  for 
stability  which  will  be  useful  for  our  purposes.  After  making  use  of 
(2.25),  (2.  27)  and  ♦ (b)  = 0 , we  can  rewrite  (6.1)  as 

b 

VU } = J C W(2  + r V2)  - W(2  + r V2)  - f(r*0(r* '-  r0O) Znrdr  (6.  5) 

a 

for  any  i/  in  X • It  follows  that  a sufficient  condition  for  the  stability 
of  0(r)  is  that  for  every  f in  x » ♦ ¥ 0 , 

W(2  + r V2)  - W(2  + r V2)  - ffr^r*'-  r0O  ^ 0 (6.6) 

at  each  r in  (a,  b)  where  the  left  hand  side  exists,  and 

W(2  + rV2)  - W(2  + r20'2)  - f(r0/)(r*/-r0/)>O  (6.  7) 

at  each  r in  some  sub-interval  of  (a,  b)  where  the  left  hand  side  exists. 
On  the  other  hand,  if  at  each  r in  (a,  b)  where  k(r)  = r0'(r)  exists, 
we  have 

W(2+k2)-  W(2  + k2(r))-f(k(r))(n  -k(r))>0  (6.8) 

for  all  numbers  K^k(r)  , it  follows  that  (6.6)  and  (6.7)  hold.  Equation 
(6.8)  is  thus  a sufficient  condition  for  the  stability  of  the  solution  0(r). 

Now  consider  the  case  in  which  both  the  inner  and  outer  surfaces 
of  the  tube  are  held  fixed  during  a virtual  displacement.  The  potential 
energy  functional  ) is  now  defined  by 
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b 

V{ ♦ } = J{W(2  + r2*'2)  - W(2  + r20/2)]2nrdr 
a 


(6.9) 


for  all  functions  i)  in  some  set  X • In  this  case  we  take  x to  be  the 
subset  of  the  previous  set  of  admissible  virtual  twists  which  conforms 
to  t(a)  = 0(a)  = 0Q  . Stability  is  defined  as  before.  By  virtue  of  (2.25) 
and  (2.27)  we  can  again  write  V{*}  in  the  form  given  by  (6.5),  whence 
(6.8)  continues  to  be  a sufficient  condition  for  stability. 


6.2  Consequences  of  Stability 

Following  Ericksen  [ 12],  we  first  make  note  of  a geometric 
property  of  the  response  curve  in  shear.  Recall  the  functions  Aj(t) 
and  A2(t)  defined  by  (5.23)  and  (5.24),  representing  the  areas  of  the 
loops  formed  by  drawing  a line  in  Fig.  1 parallel  to  the  k-axis  at  a 
distance  t , Tmin^  f * , above  it.  Recall  also  that 


Al(Tc)  = A2(tc)  ' 


A,<t»Aj(t)  for  TminST<T(;.  A,(T)<A2(T)  for  %<TSTm.x  • <*•  10> 

Keeping  this  in  mind,  one  observes  the  following  properties  of  f(k) 
upon  examining  its  graph  (Fig.  1).  If  we  set 


k3-FI«Tc)  ' k4=F3(Tc)  • 

then  we  may  observe  first  that 

(i)  if  k is  any  number  such  that  either 

| k | < k3  or  | k | > k4  , 

then 


(6. 11) 


(6.  12) 
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K 

Jf(5)d5>f(k)(x  - k)  for  all  H^k  . (6.13) 

k 

Equation  (6.  13)  is  a statement  of  the  geometric  observation  that,  pro- 
vided (6.  12)  holds,  the  area  under  the  response  curve  from  k to  k , 
for  any  K^k  is  greater  than  the  area  of  the  rectangle  of  the  same  width 
and  of  height  f(k).  By  virtue  of  (2.24),  we  can  write  (6.  13)  as 

W(2+k2)- W(2  + k2)-f(k)(K -k)>0  for  all  n^k  . (6.14) 

Next,  it  may  be  noted  that 

(ii)  if  k is  any  number  such  that 

k3<  |k|<k4  , (6.  15) 

then  there  exists  some  sub-interval  J of  (-oo,  oo)  such  that 

K 

Jf(5)d5<f(k)(n  - k)  for  h in  «9  , (6.16) 

k 

whence,  by  (2.  24),  we  have 

W(2+k2)  - W(2  + k2)- f(k)(K  - k)<0  for  x in  J . (6.17) 

Alternatively,  (i)  and  (ii)  can  be  established  analytically. 

We  now  conclude,  by  virtue  of  (6.8),  (6.  12)  and  (6.  14),  that 
any  equilibrium  solution  0(r)  for  which  either 

|k(r)|  = |r0'(r)|<k3  or  |k(r)  | = | r$'(r)  | >k4  , (6.18) 

at  each  r in  (a,b)  where  0'  exists,  is  stable.  It  is  a trivial  exercise 
to  examine  all  the  available  equilibrium  solutions  0(r)  — given  by  (4.  17)  — 
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and  to  determine  those  that  conform  to  (6.  18).  One  finds  that  only  the 
following  do: 

b 

(i)  Smooth  Solution  1 with  I 0q  I s J ^ Fj(a^Tc/5^)d5  , 

a 


b 

(ii)  Smooth  Solution  3 with  I I a J j F^(b^Tc/52)d5  , 

a 

(iii)  Weak  Solution  (3,  1)  with  |t|  = 2ir"r^Tc  , i.e.  the  solution  {3,  1)  with 

the  torque  given  by  | | =03j(.\/!  T|  / 2ittc  , | T | ) and  with  shock 

radius  7 = ^|  t|  /2tttc  . 

These  solutions  therefore  are  stable.  On  Fig.  9,  these  refer  to  the 
solutions  associated  with  points  on  the  curves  OX,  YO'  and  XY  respec- 
tively. 

We  will  now  show  that  all  the  other  solutions  are  unstable.  We  do 
this  by  exhibiting  particular  admissible  functions  *(r)  which  render 
Vf*}  negative.  It  is  readily  shown  that  these  remaining  solutions  — 
i.e.  solutions  (4.  17)  which  do  not  conform  to  (6.18)  — all  have 


k3<  |r*'(r)|<k4  , (6.19) 

on  some  sub-interval  of  [a,b].  In  each  case  \)(r)  is  chosen  to  take 
advantage  of  (6.  15),  (6.  17)  and  (6.  19).  We  first  consider  the  case  of 
dead  loading. 

Consider  for  example  Solution  1 with 


b 


b 


(6.20) 
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Recall  that  in  this  case 


(6.21) 


It  follows  from  (6.20),  (6.21)  because  of  the  monotonicity  of  the  function 


F j that 


2wa2T  <|T|s2ira^T 

c max 


(6.22) 


By  virtue  of  this,  there  is  a number  s , a<s<b,  such  that 


T|>2Trr  t for  a£r<s  . 
c 


(6.23) 


It  is  readily  shown  that  k4>  | r$'(r)  | > k3  on  (a,  s).  Let  6 beany 
number,  0<6s  s-a.  We  now  choose  the  function  t(r)  in  x such  that 


r 


Fj( — for  a<r<a+6  , 
2irr 


r*'(rH 


-Fj^ — for  a + 6<r<b  . 
V.  2wr 


(6.24) 


Observe  from  (6.2),  (6.21)  and  (6.24)  that 

a+6 


1 J 2irr  • 


2irr 


,r f i 

f_T> 

1 F 1 

f_L> 

,t 

W J 

l-Fjl 
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Therefore,  given  any  number  «>0  , we  can  choose  6 sufficiently 
small,  0<6<s-a,  such  that  ||t-d||<«. 

On  using  (6.21)  and  (6.24)  in  (6.5)  we  find  by  virtue  of  (2.24), 
(2.28),  (2.34),  (2.40)  and  (2.41)  that 

a+6  F3(r,) 

Vf  ♦ 3 = J 2nrJ  ff(5)  - n}d5dr  , r|  = -L^  . (6.26) 

a Fj(ti)  2irr 

This  can  in  turn  be  written  as 

a+6 

Vff}  = -J  A(-l-^2irrdr  (6.27) 

a 2irr 

because  of  (5.23),  (5.24)  and  (5.39).  Since  (6.22)  and  (6.23)  imply  that 

Tmax*~i>Tc  £or  aSr<a+6  . (6.28) 

2wr 

it  now  follows  from  (5.43),  (6.  27)  and  (6.  28)  that 

VU)<0  . (6.29) 

Therefore  Solution  1 with  (6.20)  in  effect  is  unstable.  On  Fig.  9 these 

solutions  are  associated  with  points  on  the  curve  XP  (excluding  X). 

The  instability  of  the  other  solutions  may  be  established  in  an 
entirely  analogous  manner,  i.e.  taking  advantage  of  (6.  16),  (6.  19) 
to  choose  t(r)  arbitrarily  close  to  *(r)  such  that  V{*}<0. 

Instability  in  the  case  when  both  boundaries  are  held  fixed  may 
be  established  in  a similar  manner,  taking  care  now  to  satisfy  the 
boundary  condition  t(a)  = d(a)  = 0Q  . For  example,  consider  Solution  1 
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with  (6.20)  in  effect.  Let  c(6)  be  the  function  defined  implicitly  by 


)d!+!  iFi(7rl)d5 


c(6) 


(6.30) 


By  virtue  of  (2.38),  (6.21)  and  the  implicit  function  theorem,  one  can 
show  that  (6.30)  does  in  fact  define  a function  c(6)  which  is  twice  * 
continuously  differentiable  in  a neighborhood  of  6 = 0 , and  that 


c(0)  = a , c/(0)>0  . 


(6.31) 


Thus  c(6)>  a for  sufficiently  small  positive  6 . 

Now  consider  the  virtual  twist  t(r)  defined  by 


r 


vJt^CttiV5  on  [a'c(6)]  • 

a 2ir5 


K. 


5 1V2ir5 


on  [c(6),b]  , 


(6.32) 


for  a sufficiently  small  6>0  . Note  that  f is  in  x by  virtue  of  (6.30). 
Observe  from  (6.2),  (6.21)  and  (6.32)  that 


= 2irarF.(— - F.(— c/(o)6+o(6)  as  6-0  (6.33) 

v 3X2tra4  1X2wa 

Twice  continuously  differentiable  when  T<  2^s^TmjLX  • The  argument 
presented  here  can  be  readily  modified  in  the  case  T = 2iri^Tm^  . 
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where  we  have  also  used  (6.31)  in  the  second  equation.  Therefore, 
given  any  number  s>0  we  can  choose  6>0  sufficiently  small  such 
that  1 1 f - <t>  1 1 < c . 

If  we  now  set 


V1(6)  = j [W(2  + rS/Z)-W(2  + r20/2)-f(r0/)(r*/-r0O}2Trrdr  , 


V2(6)  = JfW(2  + rV2)  - W(2  + r20/2)  -ffr^fr*'-  rtfO^rdr  , 

c(  6) 


(6.34) 


with  t given  by  (6.32),  we  may  write  (6.5)  as 


Vf*}  = V1(6)  + V2(6)  . 


(6.35) 


We  find  from  (6.21),  (6.31),  (6.32)  and  (6.34)  that 


V2(0)  = 0 , V/2(0)  = 0 , 


(6.36) 


because  of  (2.24)  and  (2.28).  Likewise  we  find 


V1(0)  = 0 , ^'(0)  = - a(-I^V(0)  , 


(6.37) 


where  we  have  also  used  (2.34),  (2.40),  (2.41),  (5.23),  (5.24),  (5.39) 
and  (6.31).  Note  because  of  (6.20),  (6.21)  and  the  monotonicity  of 
that  Tc*<  | X| /Zrra2^  Tmaoc  * whence  bY  (5*  43),  (6.  31)  and  (6. 37)  we  have 


v1,(0)  = -a(-^V(0)<o  . 


(6.38) 


On  using  (6.36),  (6.37)  in  (6.35)  we  find 
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Vf ♦ } = - a(-*-2-^N)c/(0)6  + o(6)  as  6-*0  , (6.39) 

so  that  by  (6.38)  V{t}<0  for  sufficiently  small  6>0  . This  establishes 
the  instability  of  Solution  1 with  (6.20)  in  effect  in  this  case.  The  instability 
of  the  other  solutions  may  be  likewise  established.  TTiis  completes  our 
instability  analysis. 

The  preceding  results  lead  to  the  conclusion  that  the  only  stable 
solutions  1 are  the  ones  given  by  (i),  (ii),  (iii)  following  Equation  (6.  18).  Recall 
that  on  the  torque-twist  diagram.  Fig.  9,  these  are  the  solutions  associated 
with  the  curve  OXYO*.  We  therefore  have  that  there  is  a unique  stable 
solution  0(r)  to  the  boundary  value  problem  in  its  weak  formulation  cor- 
responding to  every  value  of  the  applied  twist  , i.e.  there  is  a unique 
solution  0{r)  to  (4.  8)  - (4. 10)  which  conforms  to  (6. 2),  (6.3),  (6.5).  Note 
that  at  every  value  of  , the  displacement  equations  of  equilibrium  are 
elliptic  on  II  (r/7)  at  this  unique  solution. 

We  now  refer  to  a remark  made  in  Section  4.  1 that  a configuration 
involving  more  than  one  elastostatic  shock  is  unstable.  In  the  case  of 
a solution  with  a single  shock  we  showed  instability  whenever  (6.  19)  held. 

Clearly,  it  is  (6.  19)  and  not  the  number  of  shocks  that  is  important  in  that 
2 

argument.  It  is  readily  established  that  (6.  19)  holds  for  every  weak 
solution  involving  more  than  one  elastostatic  shock.  This  is  most  easily 
seen  from  a visualization  of  such  a solution  in  the  manner  explained  in 

1 These  solutions  exist  irrespective  of  the  geometric  and  constitutive 
details,  i.e.  even  in  the  cases  when  (3.23)  does  not  hold  these  are  the 
only  stable  solutions. 

2 

As  remarked  previously,  the  importance  of  (6.  19)  for  instability  is 
related  to  the  property  (6.  15),  (6.  17). 


association  with  Fig.  7.  This  leads  to  the  instability  of  such  a solution. 

Now  consider  a quasi- static  loading  of  the  body.  If  the  loading 
is  performed  in  a manner  which  involves  disturbances  of  the  type  we 
have  allowed  for,  one  would  expect  that  at  each  instant,  the  tube  seeks 
out  the  configuration  which  is  stable  against  such  disturbances.  On 
increasing  the  applied  twist  we  might  then  expect  the  loading  path  to 
follow  the  curve  OXYO' in  Fig.  9.  We  would  first  observe  a smooth 
configuration  of  the  tube.  An  elastostatic  shock  would  then  emerge  at  the 
inner  boundary  and  gradually  move  outwards,  dissappearing  upon  reach- 
ing the  outer  boundary  and  giving  way  to  a smooth  configuration.  On 
decreasing  the  applied  twist,  we  would  observe  this  process  in  reverse. 
Note  from  Section  5.  2 that  this  loading  path  conforms  to  the  dissipation 
inequality,  even  though  we  did  not  demand  — here  — that  it  should  be  so. 

In  fact  this  is  the  dissipation -free  path  referred  to  previously. 
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